2D PCA PLOT

rm(list=ls())
library(rgl)
library(dplyr)
library(tidyr)
library(plotly)

axis = list(showline=FALSE,
            zeroline=FALSE,
            gridcolor='#ffff',
            ticklen=4,
            titlefont=list(size=13))


eigenvec <- read.table("/home/gill/Project_2022/PCA_analysis/pca1.eigenvec")
eigenvec <- eigenvec[,2:22]
eigenvec <- separate(eigenvec,V2,c("Location","ID"),sep="_")
eigenvec$ID <- as.factor(eigenvec$ID)
fig <- eigenvec %>% plot_ly() %>% add_trace(type='splom',
                                            dimensions=list(
                                              list(label='PC1', values=~V3),
                                              list(label='PC2', values=~V4),
                                              list(label='PC3', values=~V5),
                                              list(label='PC4', values=~V6)
                                            ),color= ~ID,colors = c('#636EFA','#EF553B','#00CC96',"purple"),
                                            marker=list(size=7,line = list(
                                              width = 1,
                                              color = 'rgb(230,230,230,230)'))
                                            )
fig <-  fig %>% style(diagonal = list(visible = FALSE)) %>% layout(
  hovermode='closest',
  dragmode= 'select',
  plot_bgcolor='rgba(240,240,240, 0.95)',
  xaxis=list(domain=NULL, showline=F, zeroline=F, gridcolor='#ffff', ticklen=4),
  yaxis=list(domain=NULL, showline=F, zeroline=F, gridcolor='#ffff', ticklen=4),
  xaxis2=axis,
  xaxis3=axis,
  xaxis4=axis,
  yaxis2=axis,
  yaxis3=axis,
  yaxis4=axis)
fig